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ABSTRACT 

We report the results from a 2011 Suzaku observation of the nearby low-ionization 
BAL quasar/ULIRG Markarian 231. These data reveal that the X-ray spectrum has 
undergone a large variation from the 2001 XMM-Newton and BeppoSAX observations. 
We interpret this finding according to a scenario whereby the X-ray continuum source 
is obscured by a two-component partial-covering absorber with Nh ^10 22 and ^10 24 
cm~ 2 , respectively. The observed spectral change is mostly explained by a progressive 
appearance of the primary continuum at < 10 kcV due to the decrease of the covering 
fraction of the denser absorption component. The properties of the X-ray obscuration 
in Mrk 231 match well with those of the X-ray shielding gas predicted by the theoretical 
models for an efficient radiatively-driven acceleration of the BAL wind. In particular, 
the X-ray absorber might be located at the extreme base of the outflow. We measure 
a 2-10 keV luminosity of L2-10 = 3.3 x 10 43 erg s _1 for the 2011 data set, i.e. an 
increase of 30% with respect to the 2001 value. 
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1 INTRODUCTION 

According to most of the models of merger-driven 
quasar /massive galaxy co-evolution, quasars should undergo 
an heavily obscured phase associated with the presence of 
massive outflows capable of sweeping the bulk of the gas 
reservoir away from the surrounding host galaxy (e.g. Silk 
& Rees 1998; Fabian 1999; Cattaneo et al. 2009; Zubovas 
& King 2012). Mrk 231 is regarded as one of the most 
promising candidates of a quasar that gradually emerges 
from the obscured phase prior to becoming X-ray unob- 
scured/optically bright, while the AGN feedback onto the 
host galaxy ISM is active. Mrk 231 is indeed unique for many 
reasons (e.g. Lipari et al. 2009). Firstly, it is the nearest 
quasar (z = 0.042) and the most luminous (Lm=Ls-ioooiim 
= 3.6 x 10 12 Lq) ULIRG in the local Universe, with a sig- 
nificant fraction (>40%) of its bolometric luminosity due to 
starburst activity (SFR ~ 200 Mq yr -1 ), with typical signa- 
tures of a late-state merger. Furthermore, Mrk 231 belongs 
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to the rare subclass of low-ionization broad absorption line 
(LoBAL) quasars, which exhibit low-ionization (i.e. Mgll 
and A1III) BALs and very weak [OIII] emission in their op- 
tical spectra. Their unusual properties were explained by 
Boroson & Meyers (1992) in terms of absorption with higher 
column density and larger covering fraction than in other 
BAL quasars, which blocks the ionizing radiation from the 
nucleus that normally illuminates the narrow- line region. Su- 
pergiant expanding bubbles of ionized gas on kpc-scale have 
been also detected (see Lipari et al. 2009). 

The remarkable discovery of spatially-extended, broad 
wings in the CO(J=1-0) emission line profile with inter- 
ferometric millimeter observations of Mrk 231 (Feruglio et 
al. 2010) revealed the presence of a giant molecular out- 
flow with an estimated mass outflow rate of ~700 Mq yr -1 , 
which largely exceeds the ongoing SFR in the host galaxy. 
This finding provided one of the very first direct and com- 
pelling evidence of quasar-driven feedback removing large 
amounts of molecular gas, i.e. the material of star forma- 
tion, from the central regions of the galaxy (then confirmed 
by Rupke & Veilleux 2011, Aalto et al. 2012 and Cicone et 
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al. 2012). All these results support the hypothesis that in 
Mrk 231 both the nucleus and the host galaxy are under- 
going a key phase, still poorly investigated, in the evolution 
from an ULIRG/gas-rich system to a standard, X-ray- and 
optically-bright quasar in a passive galaxy. 

This motivates to shed further light on the properties of 
the central engine in Mrk 231 via penetrating X-ray obser- 
vations. However, despite its peculiarity and distinctiveness, 
Mrk 231 has not been targeted by an extensive, well-defined 
program of X-ray observations so far. The ASCA data taken 
in 1999 suggested the presence of multiple reprocessed emis- 
sion components as possible explanation for the very com- 
plex spectrum (Maloney & Reynolds 2000). Then, four short 
(40 ks each) exposures have been performed by Chandra 
(Gallagher et al. 2005), while XMM-Newton observed this 
source for only ~20 ks in 2001. Braito et al. (2004, B04 
hereafter) published these data along with those obtained 
by BeppoSAX a few months later, revealing the presence of 
a 3cr hard X-ray (> 15 keV) excess in the Phoswich Detector 
System (PDS) spectrum, suggesting that the X-ray primary 
continuum is heavily obscured by a Compton-thick (CT) ab- 
sorber (i.e. Nh Si 1.6 x 10 24 cm -2 ). This matches well with 
X-ray studies of LoBAL quasars that typically find highly 
absorbed spectra (Morabito et al. 2011; Ballo et al. 2011). 
B04 derive a quasar-like intrinsic luminosity of L2-10 ~ 10 44 
erg s _1 for the AGN in Mrk 231. However, it is important 
to stress the large uncertainties affecting these data and, 
therefore, the inferred spectral parameters above 10 keV. 

Here, we present the first 0.7-30 keV observation of Mrk 
231 obtained with Suzaku. The unique broad-band coverage 
of Suzaku allows the simultaneous study of the complex, 
absorption dominated-spectrum below 10 keV, and the very 
hard X-ray featureless continuum of the AGN, free from con- 
tamination from absorption. The Suzaku data confirm the 
BeppoSAX/PDS detection of a high-energy excess emission 
above the extrapolation of the 2-10 keV spectrum. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 231 has been observed by Suzaku (Mitsuda et al. 2007) 
on 2011 April 27 (Obs. ID: 706037010) for 200 ks. During 
the exposure, X-ray Imaging Spectrometer (XIS) CCDs (e.g. 
Koyama et al. 2007) and Hard X-ray Detector (HXD; e.g. 
Takahashi et al. 2007) have been operated in their normal 
mode, with both front- and back-illuminated XIS chips in 
the 3x3 and 5x5 editing modes. Owing to an instrumental 
failure, XIS2 is unavailable for observations performed after 
November 2006. To maximize the HXD effective area, the 
target has been observed at the HXD nominal aim point. 
Data from XIS and HXD have been reduced with FTOOLS 
6.12 and Suzaku software Version 19, using the latest calibra- 
tion files available (2012 February 11 release) and following 
the prescriptions described in the Suzaku Data Reduction 
GuidfS 

For the XIS detectors (i.e. XIS0, 1 and 3) source counts 
are extracted from a circular source-free region with a radius 
of 146 arcsec centered on the peak of the X-ray emission. 
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Background counts are estimated from a similar circular re- 
gion located close to the source. The XIS response and an- 
cillary response files have been produced using the FTOOLS 
tasks xisrmfgen and xissimarfgen, respectively. The 3x3 and 
5x5 spectra and matrices have then been combined together 
using addascaspec. All spectra are binned to a minimum of 
30 counts per bin to allow use of the \ 2 statistic. During the 
Suzaku observations Mrk 231 did not show any variation in 
X-ray flux and the spectral shape. 

The HXD/PIN cleaned event file and the corre- 
sponding "tuned" non- X-ray background (NXB) event file 
(ae706037010_hxd_pinbgd.evt) have been used as input to 
the hxdpinxbpi task. The cosmic background spectrum is 
created with a model provided by the Suzaku team based 
on the formula by Gruber et al. (1999) using a flat response 
(ae_hxd_pinflatel0_20101013.rsp) and then combined with 
the NXB spectrum extracted from the NXB event file. The 
resulting combination is used for the background subtrac- 
tion. We consider only data in the 15-30 keV energy range, 
where the source contribution is ~6% of the total signal of 
the PIN. The source contribution drops down to 3% in the 
30-50 keV energy bin. This choice ensures that the source 
signal is well above the current systematic uncertainty on 
the PIN background, i.e. ~ 3% for observation with exposure 
longer than 10 ks (Fukazawa et al. 2009). We use the appro- 
priate response file (ae_hxd_pinhxnomel0_20101013.rsp) for 
HXD-nominal pointing for the spectral analysis. Here, we 
present data from PIN only, since no significant signal was 
detected from GSO. Relative normalization among XISO-3 
and PIN are set as reported in Maeda et al. (2008) for an 
observation performed at the HXD nominal position. 

XMM-Newton and BeppoSAX observations of Mrk 231 
performed in 2001 are also used here. XMM-Newton data 
have been reduced with SAS 11.0 (Gabriel et al. 2004). 
The pn source and background spectrum and the matri- 
ces are created using standard procedures (e.g., Piconcelli 
et al. 2004). The BeppoSAX/PDS spectrum and the corre- 
sponding response matrix have been downloaded from the 
BeppoSAX public archive. 

Spectra are analyzed with the XSPEC vl2.7 package 
(Arnaud 1996). All fitting models presented here include ab- 
sorption caused by the line-of-sight Galactic column density 
of N H = 1.3 x 10 20 cm -2 (Dickey & Lockman 1990). 



3 SPECTRAL VARIABILITY AND 
INTERPRETATION 

B04 found a best-fit model (indicated as model A in their 
paper) to the joint 2001 XMM-Newton and BeppoSAX spec- 
trum of Mrk 231 consisting of a combination of a power law 
modified by a fully-covering CT absorber, and a mix of re- 
processed AGN emission absorbed by an additional screen 
with Nh ~ 10 22 cm -2 , and multi-temperature thermal emis- 
sion due to the starburst activity, which accounts for the 
very complex 0.5-10 keV spectral shape. 

Fig. la shows the application of this fitting model to 
the Suzaku data taken in 2011. This comparison implies a 
very significant X-ray spectral variation in a 10 yr time span, 
accompanied with a moderate increase of the flux over the 
Suzaku bandpass. In particular, the XIS spectral data in the 
2-10 keV band show the largest differences with respect to 
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Energy (keV) Energy (keV) 

Figure 1. (a)-Left: the unfolded Suzaku XISO (black points) + PIN (red points) spectrum of Mrk 231 collected in 2011 compared with 
the best-fit model to the XMM-Newton + BeppoSAX data taken in 2001. (h)-Right: the Suzaku spectrum fitted with the best-fit 'double 
partial-covering absorber' model (XISO: black points; XIS1: red points; XIS3: green points; PIN: blue points). The lower panel shows the 
deviations of the observed data from the model in units of standard deviation. 



the XMM-Newton ones. Moreover, Fig. la indicates that the 
normalization of the Fe Ka emission line at ~6 keV remains 
constant between the two epochs (see Sect. 4), suggesting 
no variation in both spectral features, i.e. fluorescence iron 
line and continuum, associated with the reflection compo- 
nent (George & Fabian 1991; Ghisellini et al. 1994), which 
is believed to originate from a distant reprocessor (typically 
the visible inner wall of the obscuring torus itself). The soft 
X-ray band in model A is entirely due to starburst emission. 
As such, no flux variability has to be expected in the soft X- 
ray band. The 2-10 keV band is instead dominated by X-ray 
reflection likely on 0.1 pc scales. Variability in response 
to continuum changes could in principle be observed on a 10 
yr timescale, but this is at odds with the apparent constancy 
of the associated Fe line intensity. Hence, the only plausible 
explanation for the soft and hard X-ray band flux variability 
is that part of the X-ray continuum is visible below 10 keV. 

We then replace the CT absorber of model A with a 
partial covering one which allows some continuum leakage 
below 10 keV. According to this scenario, the primary con- 
tinuum power law is seen through two distinct, nested lay- 
ers of absorption. The underlying idea is that the observed 
variation in the Suzaku spectrum with respect to the XMM- 
Newton one can be explained in terms of changes in the 
physical and geometrical properties of these absorbers. In 
our analysis, we use the zxipcf model in XSPEC (Reeves 
et al. 2008), which describes partial covering absorption by 
partially ionized material at the redshift of the source. This 
allows us to explore a large variety of possible ionization 
states and covering fraction of the obscuring medium. To 
account for the reflection component, we use the pexmon 
model (Nandra et al. 2007) which includes both continuum 
reflection (Magdziarz & Zdziarski 1995) and the associated 
Fe K line emission (George & Fabian 1991), with the metal 
abundances of the reflector fixed to the solar value and its 
inclination angle fixed to 60° . The photon index of the X-ray 
continuum and reflection is fixed to T = 1.9, i.e. the typical 
value measured for quasars (Piconcelli et al. 2005). Similarly 
to B04, we also include in the fit two mekal thermal compo- 
nents with solar abundance. This model yields an excellent 
description of the 2001 XMM-Newton+ BeppoSAX data, i.e. 
X?(dof) = 0.96(76) (Table [l]). The spectral complexity and 



the quality of the data above 10 keV do not allow to put 
tight constraints on both the column density and the ion- 
ization parameter £ (expressed in erg cm s ; see Reeves et 
al. 2008) of the this absorber. Therefore, we fix the parame- 
ters to their best-fit values deduced by the fit, i.e., Nh = 1.6 
x 10 24 cm -2 and log £ = -2. The covering fraction of this 
CT absorber is ~0.85. The covering fraction of the second 
partially-covering absorber is ~ 0.92 (while fixing C/ to the 
unity worsens the fit by Ax 2 = 26 for 1 dof). The best-fit 
values for the column density and the ionization parameter 
of this obscuring component are Nh = 5.9 x 10 22 cm -2 and 
log £ = -0.4, respectively. 

In this fit, we assume that the reflection component is 
not obscured by the two-component absorber and its nor- 
malization Ar is equal to that of the primary continuum 
Apl- However, a fit with the reflection component obscured 
only by the outer/Compton-thin absorber produces similar 
results. Nonetheless, we prefer a solution in which the reflec- 
tion component is left outside the two-component partial- 
covering absorber as the constant Fe K line flux suggests 
that the reflection originates far away the continuum source, 
while the partial-covering and warm nature of the Compton- 
thin absorption is more consistent with an inner, nuclear 
location rather than a pc-scale one. The value of the re- 
flection fraction is R ~ 0.2 (defined as Q/27T, where Q is 
the solid angle subtended by the reflector for isotropic inci- 
dent emission). The temperature of the two mekal compo- 
nents accounting for the thermal diffuse emission is fcT = 
0.4±0.1 and fcT = 1.0±0.2 keV, respectively, i.e. consistent 
with the values derived by spatially-resolved Chandra spec- 
troscopy (Gallagher et al. 2002). Finally, this model implies 
a 2-10(15-30) keV observed luminosity of 2.5(11.9) x 10 42 
erg s -1 . 



4 BROAD-BAND SUZAKU SPECTROSCOPY 

As the next step, we apply the double partial-covering ab- 
sorber (DPCA hereafter) model that has been found to suc- 
cessfully fit the 2001 XMM-Newton+ BeppoSAX data to the 
broad-band Suzaku spectrum taken in 2011. Our goal is to 
test if both data sets can be described by the same model, 
and the observed variability can be ascribed to variations 
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Table 1. Best-fitting spectral parameters for broad-band X-ray spectra of Mrk 231 measured assuming the DPCA model (Sect. [3] and 
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(2) 


N H /log i 
(3) 


Of 
(4) 


N H /log5 
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(7) 


Flux 
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XMM(6/2001) + 
SAX(12/200X) 
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160 t /-2 t 


85+ 05 
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5 Q+ 1 ' 2 / n 4 +0 - 7 
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1.1/6.4/30.2 


0.96(76) 


Suzaku{A/2011) 
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°- a -0.6/ U - a -0.4 


0.95±0.01 
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1.7/10.1/37.1 
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Notes: (1) X-ray telescope (epoch); (2) Intensity of the power-law primary continuum (in units of 10 — 3 photons keV - 1 cm -2 s — 1 ); 
(3) Column density (in units of 10 22 cm -2 ) and ionization parameter of the CT absorber; (4) Covering fraction of the CT absorber (5) 
Column density (in units of 10 22 cm -2 ) and ionization parameter of the Compton-thin absorber; (6) Covering fraction of the Compton- 
thin absorber; (7) Reflection fraction; (8) Flux in the 0.5-2/2-10/15—30 keV band (in units of 10 — 13 erg cm -2 s _1 ). t The parameter 
has been fixed to the best-fitting value deduced by the fit. 



in the line-of-sight absorbers. The normalization of Fe Kq 
emission line measured from the Suzaku spectrum is Ap e 
= 9±6 x 10 -7 photons cm -2 s _1 . Such a value is fully 
consistent with the XMM-Newton one (Af e = 1.5±1.1 x 
10~ 6 photons cm -2 s _1 ) and those inferred from Chandra 
observations (e.g. Gallagher et al. 2005). The Fe Kq line 
and, hence, the associated reflection continuum therefore re- 
mains approximately unchanged over these epochs (see Fig. 
la). Accordingly, the value of reflection fraction R has been 
forced to vary within the 90% confidence level interval de- 
rived for these parameter from the 2001 data, keeping the 
same normalization (Ar = 2 x 10~ 3 photons keV -1 cm" 2 
s _1 ). We then fix R = 0.08, i.e., the best-fit value deduced 
by the fit. The values of the Nh and £ of the CT absorber 
are constrained as for the XMM-Newton+BeppoSAX spec- 
trum. The result of the DPCA fit to the 0.7-30 keV Suzaku 
data is shown in Fig. lb. This model provides a very good 
fit with x^(dof) = 1.09(285). 



As expected on the basis of Fig. QJ,, we measure an in- 
crease of the flux in the 0.5-2, 2-10 and 15-30 keV bands (e.g. 
Table[l} from 2001 to 2011. It mostly results from the com- 
bination of a larger normalization of the power law of the 
primary continuum and a smaller covering fraction of the 
CT absorber in the Suzaku spectrum (although, given their 
large uncertainties, both parameters are however consistent 
at 90% with those inferred in Sect. [3]). The best-fitting tem- 
peratures (AT = 0.23+0^ and AT = 1.26±g;|| keV) and the 
normalizations of the two mekal soft X-ray emission com- 
ponents are in good agreement with those derived from the 
XMM-Newton spectrum, as expected for spectral compo- 
nent associated with starburst emission and/or photoionized 
gas emission. Once model DPCA is assumed, we estimate a 
2-10(0.5-30) keV flux of ~1(5.7) x 10" 12 erg cm -2 s" 1 , cor- 
responding to an observed luminosity of ~4(22) x 10 42 erg 
s" 1 . 



Finally, replacing the warm absorbers in the DPCA 
model with cold ones yields a worse fit to the Suzaku data, 
with an associated xt(doi) = 1.28(285), and large positive 
residuals around 5 keV where it predicts less spectral cur- 
vature and a continuum less intense than observed. 



5 DISCUSSION 

The analysis of the Suzaku data presented here has revealed 
a large variation of the X-ray spectral shape from the last X- 
ray observations of Mrk 231 performed with XMM-Newton 
and BeppoSAX in 2001 (Fig. UK). We have successfully fitted 
the broad-band Suzaku spectrum with a DPCA model. This 
model also provides an excellent fit to the 2001 data. In this 
scenario, two absorbers lie along our line of sight (LoS) to 
the nucleus of Mrk 231: one with a CT column density, and 
another one with Nh ~ a few of 10 22 cm -2 . Both absorbers 
are slightly ionized and have a covering fraction less than 1. 
This allows to detect the leakage below 10 keV of the direct 
emission from the X-ray source. Accordingly, the spectral 
variation shown in Fig. [TK can be basically explained by a 
progressive appearance of the primary continuum at energies 
<10 keV due to a decrease of the covering fraction of the CT 
absorber, suggesting that this absorption component is the 
innermost, more variable one. This interpretation appears to 
be consistent both with (i) the spectral properties inferred 
by previous X-ray observations, as well as (ii) the LoBAL 
nature of the quasar in Mrk 231. 

Maloney & Reynolds (2000) indeed proposed the ex- 
istence of two layers of obscuration after the analysis of a 
moderately long ASCA exposure. However, they assumed 
that no fraction of the X-ray primary continuum is visible 
in the ASCA range, i.e. the quasar central engine is com- 
pletely buried. If, as we suggest, the C/ variation occurred 
between 2001 and 2011, it is indeed possible that ASCA 
caught the source at an epoch when the C/ of the CT ab- 
sorber was close to unity. However, it is important to bear 
in mind that the limited bandpass of ASCA did not allow 
to put any meaningful constraint on the CT obscured X- 
ray continuum. The variability in the hard X-ray flux de- 
tected during a 40 ks Chandra observation led Gallagher 
et al. (2002) to question the reflection-dominated hypoth- 
esis for the X-ray spectrum of this quasar, although such 
variability was not detected in any of the subsequent three 
similar Chandra exposures performed in 2003. They sug- 
gested the existence of a third spectral component, i.e. an 
absorbed (Nh ~ 3 x 10 22 cm -2 ) power law, in addition to 
the thermal plasma and the reflection continuum, which may 
arise from a compact X-ray scattering medium. Gallagher 
et al. (2005) found an acceptable fit to the combined 0.5- 
8 keV spectrum from the four Chandra observations with a 
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model consisting of a power law modified by three absorbers 
plus soft thermal emission, which may resemble our DPCA 
model in terms of overall spectral shape in this energy band, 
'weak' Fe Ka emission and, mostly, for the presence of a 
multi-component absorber. Furthermore, the consistency of 
the flux of the Fe Ka emission line over 2001-2003 found 
by Chandra and XMM-Newton observations of Mrk 231 is 
suggestive of a Compton reflection component likely to be 
associated with distant, optically thick material, and, re- 
markably, the 2011 Suzaku data lend further support to this 
scenario. The DPCA model proposed here is therefore able 
to explain in a coherent fashion the several pieces of evi- 
dence collected over the years thanks to the previous limited 
bandpass and/or less sensitive X-ray observations. The hard 
X-ray flux values measured from the Chandra observations 
are consistent each other and are slightly larger (i.e. consis- 
tent within 2a) than the XMM-Newton value (Gallagher et 
al. 2005). In the proposed DPCA scenario, this can be inter- 
preted as lack of large variations in the physical properties 
of the CT absorber. 

The presence of the highly absorbed power-law emis- 
sion at E > 15 keV above the extrapolation of the 2-10 keV 
spectrum claimed by B04 thanks to the 2001 BeppoSAX 
observation is fully confirmed by the present Suzaku data. 
X-ray spectral properties of LoBAL are very scarce due to 
their weakness in the 0.5-10 keV range (Green et al. 2001; 
Morabito et al. 2011), interpreted in terms of heavy obscu- 
ration. Even much more extraordinary is the detection of 
these objects at higher energies: as far we know, Mrk 231 is 
the only member of this class of AGN for which a spectrum 
has been collected in the band above 10 keV. 

The change in the spectral shape discovered here as well 
as the short-timescale variability in the hard-band emission 
reported by Gallagher et al. (2002) (assuming an origin in 
the multi-component absorber) seems to indicate that the 
bulk of the X-ray obscuration likely occurs in a compact re- 
gion much closer to the X-ray source than the canonical 0.1- 
1 pc-scale torus invoked in the Unified model, or extended 
structures in the complex circumnuclear environment of Mrk 
231 (e.g. Klockner et al. 2003, and references therein). Both 
ionized and neutral X-ray absorbers with large column den- 
sities are commonly observed in BAL quasars (e.g. Green 
et al. 2001; Piconcelli et al. 2005; Morabito et al. 2011; 
Giustini et al. 2011) They have been indeed associated with 
the BAL outflow (e.g. Elvis 2000), which is likely launched 
from the accretion disk surrounding the supermassive black 
hole (Young et al. 2007), and their variability explained in 
terms of outflowing clouds crossing our LoS. In particular, 
it has been proposed that this wind must be shielded from 
the X-ray continuum source to prevent the BAL absorb- 
ing clouds from becoming too highly ionized. This would 
inhibit a very efficient mechanism of gas acceleration such 
as the radiation-driven acceleration, by eliminating most of 
the UV transitions that absorb the photons necessary to ac- 
celerate the outflowing gas (Murray et al. 1995; Proga 2007, 
and references therein). This X-ray-shielding region must be 
effectively opaque to X-rays, as well as sufficiently ionized 
(and dust-free) to allow UV line opacity, and might rep- 
resent the inner section of the wind itself (and/or located 
at the extreme base of the outflow). The properties of the 
DPCA discovered in Mrk 231 match well with those pre- 
dicted for the shielding gas by the BAL outflow models, 



providing an important observation-based evidence for the 
theory of physical mechanism responsible for launching and 
accelerating AGN winds. These nuclear, highly supersonic 
outflows are believed to compress the surrounding medium 
and generate a shock, which moves outwards and sweeps the 
surrounding medium out, i.e. the AGN feedback mechanism 
(Menci et al. 2008; Zubovas & King 2012). 

The presence of a partially-covering absorber suggests 
that the typical size of the X-ray absorber along our LoS is 
comparable with that of the X-ray continuum source. The 
latter is likely of the order of ~10 Re, as estimated via oc- 
cupation or microlensing in a few AGNs (Dai et al. 2010; 
Risaliti et al. 2007; Maiolino et al. 2010). Accordingly, the 
DPCA is not the torus envisaged by the Unification model, 
which can be however present in Mrk 231, as it should lie 
(i) on larger scales and (ii) edge-on with respect to our LoS. 
The non-variable reflection emission component observed in 
the X-ray spectrum of Mrk 231 is likely associated to this ob- 
scuring/reprocessing medium located distant from the cen- 
tral engine. 

Finally, the broad bandpass of Suzaku allows us to mea- 
sure the intrinsic luminosity of the X-ray source in Mrk 231. 
We find an intrinsic luminosity of L2-10 = 3.3 x 10 43 erg 
s _1 in the 2-10 keV energy range, while the luminosity in the 
15-30 keV PIN band is 1.7 x 10 43 erg s" 1 . Using the DPCA 
model for the joint XMM-Newton + BeppoSAX spectrum, 
we calculate a L2-10 = 2.6 x 10 43 erg s _1 . It is interest- 
ing to compare this X-ray luminosity with that expected 
on the basis of well-established isotropic indicator of AGN 
luminosity. The luminosity of the mid-IR emission line [O 
IV] 25.89/im has been shown to be one of the most accurate 
proxies of L2-10 (Melendez et al. 2008). Armus et al. (2007) 
reported an upper limit for the [O IV] 25.89/im luminosity 
of i[oiv] < 2.64 x 10 41 erg s" 1 based on 2004 Spitzer data. 
Using the L[ iv] _ i2~io relation in Melendez et al. (see Ta- 
ble 3), the predicted 2-10 keV luminosity for Mrk 231 is < 
3.7 x 10 43 erg s _1 . This matches well with the L2-10 values 
deduced from both Suzaku and XMM-Newton observations. 
Similar results are also obtained by the comparison of the 
luminosity in the 14-195 keV range predicted on the basis of 
Lroiv], i- e - £14-195 < 10 44 erg s _1 , and those estimated from 
Suzaku and XMM-Newton + BeppoSAX data, i.e. L14-195 
= 7 and 5.7 x 10 43 erg s _1 , respectively. 

Future X-ray observations with broad-band X-ray capa- 
bilities (i.e., NuSTAR and ASTRO-H) as well as well-defined 
monitoring program with present-day X-ray telescopes will 
allow to provide much more detailed insights into the ki- 
netic, geometrical and physical properties of the nuclear ab- 
sorber^) in this composite ULIRG/LoBAL quasar, where 
AGN-driven outflows in the neutral, ionized and molecular 
gas capable of quenching star formation in the host galaxy 
have been found. 
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